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Abstract 



I'T I We analyze the non-equilibrium steady states (NESS) of a one dimensional harmonic chain of N 

^ atoms with alternating masses connected to heat reservoirs at unequal temperatures. We find that 

^ the temperature profile defined through the local kinetic energy T{j) = < p^- >/mj, oscillates with 

Qh period two in the bulk of the system. Depending on boundary conditions, either the heavier or 

I 

^5 the lighter particles in the bulk are hotter. We obtain exact expressions for the bulk temperature 

^ profile and steady state current in the limit N ^ oo. These depend on whether N is odd or even. 

We also study similar temperature oscillations in the NESS of systems with noise in the dynamics. 
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I. INTRODUCTION 



The study of non-equilibrium steady states (NESS) of macroscopic systems in contact 
with heat baths at different temperatures has a long history jTH2]. There are no known 
analytic solutions for interacting Hamiltonian systems, except for harmonic crystals. When 
the atoms in the crystal all have the same mass, this NESS can be obtained explicitly 
[His]. It gives a uniform "temperature", i.e (p|)/m, in the bulk of the system. It also gives 
heat currents that are independent of the size of the system corresponding to the fact that 
phonons can travel freely through the crystal. This behavior of the heat current is also true 
for harmonic systems with periodic arrangement of masses [HI E] • It was therefore surprising 
when numerical simulations of the NESS showed that the bulk temperature profile of a 
chain with alternating masses oscillates between two values and that these oscillations did 
not seem to decay on increasing the system size [Ij. 

Here we present analytical solutions for the temperature profile and current in the 
alternate-mass chain connected to Langevin type heat baths, which prove that the oscil- 
lations persist in the — t- oo limit. Only for very special choice of parameter values, can 
the oscillations be made to vanish. Surprisingly, the values of the oscillating temperature 
and of the current depends on whether is even or odd even in the asymptotic system size 
limit. An oscillating temperature profile in a thermodynamically large system is surprising 
when we think from the standard view-point of heat fiow occuring from hot to cold regions. 
However this expectation will be true only in systems exhibiting local thermal equlibrium 
where one can define a meaningful thermodynamic local temperature. This is the case for 
a system with normal diffusive heat transport, though a microscopic derivation of the con- 
ditions when this is achieved is in general difficult [2j. In the study of systems in NESS 
it is natural to define a local "temperature" from the mean local kinetic energy and this is 
what we do here — the absence of local equilibrium in the harmonic chain allows for the 
"temperature" profile to show the unexpected oscillatory feature. 

In the context of testing Fourier's law and investigating the size-dependence of the current, 
various one-dimensional models with alternate masses have been studied P-IIO). Tempera- 
ture oscillations have earlier been observed in the steady state of the alternate mass hard 
particle gas [9J and in the Fermi-Pasta-Ulam chain [lOJ but in these cases the oscillations de- 
cay with system size. The case of temperature oscillations persisting for infinite system sizes 
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is thus special to harmonic systems where heat is transmitted by non-interacting phonons. 
It is expected that introduction of phonon-phonon interactions will in general make things 
different. Here we investigate this issue by considering alternate-mass harmonic chains where 
the dynamics is stochasticlly perturbed by noise which either conserves both energy and mo- 
mentum or conserves only energy. Finally to consider the effect of dimensionality, we present 
results from simulations of two-dimensional strips of alternate mass harmonic systems. 
The plan of the paper is as follows. In Sec. (|TT| we define the precise model and present 



some of the numerical results for small finite systems. In Sec. (Ill) we present the analytic 



and numerical results in the limit N — )► oo. In Sec. (IV) we present simulation results on 



temperature profiles in harmonic chains with noisy dynamics. In Sec. (|Vj) we summarize 
our results and give a physical explanation of the results. The details of our analytical 
calculations are given in Appendix (lAl). 



II. MODEL AND NUMERICAL RESULTS FOR SMALL SYSTEM SIZES 

We consider a one-dimensional chain of particles labeled i = 1, . . . , N that are placed 
in an external harmonic potential (with spring constant k^) and which are interacting with 
each other through a nearest neighbor harmonic potential (with spring constant k). Let the 
vectors q = {qi, q2, . . . , ^iv) and p = {pi,P2, ■ ■ ■ ,Pn) denote respectively the displacement 
and momenta of the A^ particles of the chain. The Hamiltonian for the ID chain we consider 
is given by: 

^ i=N N+1 N 

H = Y^p\l2m, + \ k{q, - + ^ E ^oQ" (1) 

1=1 1=1 i=l 

= ^p.M-\p+^q.^.q , (2) 

with go = Qn+i = 0, and in the second line we have used a compact notation with M defining 
the mass matrix and $ the force-matrix. The ends of the chain are coupled to Langevin 
reservoirs at temperatures Tl and Tr. The equations of motion of the system is given by: 

ruiqi = - J2 '^i^J^lj + kii-^LQi + + Si,N[-lRqN + {^jHTRY^^m] , (3) 

j=l,N 

for i = 1,2, . . . , N , where i]L,rif> are Gaussian white noises chosen from distributions with 

averages {rjLif)) = {VR(t)) = ^^'^ correlations {riL{t)riL(t')) = {i]Rit)VRi't')) = 5{t — t') and 
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iLi'^R are dissipation constants ( Note that the derivation in [6] uses a different convention 
for the reservoir couphng. The dissipative forces on the end particles were there taken 
to be — AiPi and —XnPn and so their couphng constants are related to ours as Aimi = 
7l, XNfnN = 7_R )• We will be interested in the case where the masses rrii alternate between 
two values and rrih on successive sites. 

Corresponding to the Langevin equations in Eq. ^ it is straightforward to write the 
Fokker-Planck equation to describe the evolution of the phase space distribution fi{x,t), 
^ = ■ ■ ■ jQNyPi,''' )Pn)- Following standard methods [12j it can be shown that the 
Fokker-Planck equation is given by: 



N 
i=l 



N 



■piU + TiVdi — 

OPi 



(4) 




where the right hand side of Eq. Q describes the interaction of the end particles with the 
heat baths and Ti^n = Tl^r and 7ijv = 1l,r- Let us define the 2N x 2N matrix 



(5) 



where F is a x diagonal matrix with Fjj = 7i5ij(5ji + 5jAr) . We also define the 2N x 2N 
matrix d with elements dij = 2'y6ij(TL6i^N+i + Tji6i^2N)- It is known that the steady state 
distribution is Gaussian ^ and given by 

27V 
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2N 

where the covariance matrix b with elements bij =< XiXj > satisfies 

a.b + b.at = d . (6) 

The solution of the linear equations, Eq.(|6|, gives us all the correlations bij and hence the 
temperature profile Tj =< > /rrii = bN+i^^^i/nii , and the current, J = — 
qi)Pi/mi) = k{bi^i^i\f^i — bi^N^i)/mi. In the equal mass case the covariance matrix for sites 
can be obtained in a fairly explicit form This seems to be difficult for the alternate mass 
case. However the matrix equations can be solved numerically for small system sizes and we 
can obtain accurate results for the temperature profile and current for these system sizes. In 
Fig. ([T| we show typical temperature profiles for alternate mass chains with even and odd 




Figure 1: (Color online) Temperature profiles for (a) system with even number of sites N = 32,64, 
and with 'Jl = IR = 1-0 and (b) system with odd number of sites = 33,65, and with = 
1.5, 7i{ = 0.5. Other parameters were set to nia = 0.75, m;, = 0.25, A; = 1,Tl = 1.5, Tr = 0.5. 
The mass of the first particle is always taken to be m^. Note that in (a), the heavier particles 
are hotter, while in (b), the lighter particles are hotter. The horizontal dashed lines indicate the 
analytic predictions for N ^ oo, from Eqs. ( TTp2 ). 



number of sites for particular choices of parameter values and A^. We see oscillations in the 
temperatures of the particles in the bulk in both the even and odd cases, and the amplitude 
of the oscillations does not seem to change with system size. In the next section, we will 
obtain expressions for the current and the bulk temperatures and show that the temperature 
ocscillations persist in the N ^ oo limit. 
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III. ANALYTICAL AND NUMERICAL RESULTS IN iV ^ oo LIMIT 



To obtain analytic results in the limit — )■ oo we follow [B] and express the covariances in 
terms of integrals over frequencies. The integrands involve elements of the following Green's 
function: 

G+ = [-Mu^ + $ - iujT]-^ . (7) 
Here we are interested in the temperature and current and these are given by [6l [TT] 

/oo roo 
dujuj''\Gl{uj)? + lRTR rfwa;' |G+v(a;)|2 , z = 1, 2, . . . , iV 
-oo J — oo 

/oo 
dujuj^ |G+^Hr (8) 
-oo 



T 1 

li = -rrii 



- 1 

TT 



We rewrite the above expressions in the following form. 



Ti — li Tl + Jj Tr 



J — -^{Tl — Tr)In , 
rriN 

where the 

h = '^ r dco co' IG+HP , k = '^ r du |G+,(a;)p . (9) 

71" J-oo 

are independent of the temperatures and Tr. Now we note that for the equilibrium case 
Tf, = Tr, we must have the same temperature at all sites, z.e Tj = T, and hence deduce the 
equality /j + = 1. Using this fact and defining = T + AT/2, Tr = T — AT/2 we can 
rewrite the equation for the temperature profile in the following form: 

T, = T + (/, - 1/2) AT = Tr + hAT. (10) 

We thus only need to evaluate the integral /j, in the limit N ^ oo. 

So far our treatment has been quite general. We now focus on the alternate mass case. 
We define the first mass to be mi = nia and the next to be and so on. Thus odd sites 
have iTLi = rria and even sites, rrii = rrih. For simplicity we only consider the unpinned 
case ko = 0. It is straightforward to extend the calculations to the case ko 7^ 0. Without 
loss of generality we can choose time and energy scales so that k = 1 and rria + mi, = 1. 
We give the details of the calculation in the appendix. The main result is that Jj can be 
written as a sum of two parts, one coming from the acoustic modes of the system and one 
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from the optical modes. We note that the mode frequencies for the acoustic and optical 
bands are respectively given by: uj^_ = (l/m(jm{,)[l — 0(g)] ,a;^ = {l/mamb)[l + (f){q)], where 
0(g) = [1 — 2'ma'mf,{l — cos q)V^^ , ,0 < g < vr. The various expressions depend on whether 
is even or odd and for these two cases corresponding to superscript E, O respectively, we 
get: 

Case (1)- N = 2L, L ^ oo: 



E 



J' 



T + AT 

dq 

T + AT 

+ dq 
Jo 

AT 



dq 



{m,ul-2f + A^lulco^\q/2) 







27r0(g) |2(7l + 7^^) - (m„7^ + m67L)tc;^| (1 + '^lIr^V) 



[■mbU_ 



2f + A^lu^_cos\q/2) 



27r0(g) |2(7l + 7^) - {rria'yR + mfe7L)a;2 | (1 + 'JlIru^'^) 



dq 



W(g/2)+7M( 



2f 



27r0(g) |2(7l + 7^^) - (m„7^ + m;,7L)w^| (1 + ^iLlR^V) 

iLVfi^ 4cos^(g/2) + -i\u?_{mau?_ - 2f 1" 

27r0(g) |2(7l + 7^?) - {ma^R + mfe7L)a;2 | (1 + 7^7^20;^) 2_ 

iLlR 



dq 



sin^g 



+ 



dq 



7r0(g) |2(7l + 7^) - {nia^R + mfe7L)cu^| (1 + 7L7iJW+) 

iLlR 



sin^g 



■7r0(g) |2(7l + 7/j) - ima-iR + m67L)a;2 | _^ ^iLlR'^V) 
where the subscript o refers to odd sites and e to even sites. 
Case (2)- N = 2L + \, L ^ 00. 



T 



o 



T 



o 



o 



T + 


AT 




dq 




T + AT 

/•TT 




dq 




AT 


/•TT 

Jo 




dq 





dq- 



iLrria 



_ Acos\q/2)+^lul{m,ujl-2f 
^IL + 7ii)7r0(g) \{maUjl - 2) + -fL^ynuKmhul - 2)| 

-fLrria 4cos^(g/2) + -flu^_{mbul - 2)^ _ 1" 

2(7L + 7i?)7r0(g) liniaujl - 2) + -fL7R^UmbOO^- - 2)| 2_ 



dq- 



iLrrib 



im,cul-2y + A^j,ulcos\q/2) 



LJo 2(7l + 7/j)7r0(g) Km^w^ _ 2) + 7^7/?^^ (m;,a;^ - 2)| 

7Lmfc (m^cu^ - 2f + Ajlulcos\q/2) _ 1" 

2(7L + 7i?)7r0(g) |(maw2 _ 2) + ^^^^Jiim^J^^ - 2)| 2_ 



(ig 



iLlR 



•$,iv?q 



iLlR 



2) 



sin'^q 



(12) 

(7l + 7/?)7r0(g) I (m^w^ - 2) + 'JlIrUJ- {rribUj'i - 2) | J 
We now present some numerical data for the two cases of even N and odd for various 



(11) 



parameter sets. When 72, = 7ij, the above integrals can be carried out exactly, see Eqs. ( A13 



A17). In other cases, we evaluated the integrals numerically (using Mathematica). 

Case (1): We consider chains with even A^ and set 7l = 7_r = 7- In Fig. ^ we plot the 
temperatures on the odd (T^) and even {T^) sites, and also the current {J^ in inset) as a 
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Figure 2: (Color online) Temperatures at odd and even sites for a chain with even number of 
particles, plotted as a function of 7 = 7/, = 7/?, = 0.75, m;, = 0.25, A; = 1,Tl = 1.5, Tr = 0.5. 
We also plot separately the contributions of the acoustic and optical modes to the temperature at 
any site. The inset shows J and also the contributions of the acoustic and optical modes. 

function of the parameter 7. We also separately plot the contributions of the acoustic and 
optical modes to the temperatures and current. We note the following features: 

(i) Depending on the value of 7, either the heavier particles (those on odd sites), or the 
lighter ones are hotter. At 7 ^ 0.41, the temperatures at the odd and even sites are equal. 

(ii) The temperature of the heavier particles gets its main contribution from the acoustic 
modes while that of the lighter particles comes mostly from the optical modes. The heat 
current is mostly carried by the acoustic modes. 

Case (2): We consider chains with odd A^. In this case, 7^, = 7^? becomes a very 
special case: the masses of the end particles being equal, this condition implies symmetry 
between the left and right reservoirs, and this leads to a uniform bulk temperature equal 
to (Tl + Tr)/2. The more typical situation is when the two couplings are different and we 
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Figure 3: (Color online) Temperatures at odd and even sites for a chain with odd number of 
particles, plotted as a function of 7/j with = 1, nia = 0.75, mf, = 0.25, k = 1,Tl = 1.5, Tr = 0.5. 
We also plot separately the contributions of the acoustic and optical modes to the temperature at 
any site. The inset shows the variation of heat current with 7/j and also the contributions of the 
acoustic and optical modes. 

consider this by setting = 1 and changing 7^. In Fig. (|3| we plot the temperatures on 
the odd (Tf ) and even (T^) sites, and also the current (J*^ in inset) as a function of the 
parameter 7^. We also separately plot the contributions of the acoustic and optical modes 
to the temperatures and current. We note the following features: 

(i) Depending on the value of jn, either the heavier particles (those on odd sites), or the 
lighter ones are hotter. At a special value of 7_r = 7l, the temperatures at the odd and even 
sites are the same. They are both equal to the mean temperature T = 1. 

(ii) As for the even case, here also we see that the temperature of the heavier particles 
gets it main contribution from the acoustic modes while that of the lighter particles comes 
mostly from the optical modes. The heat current is again mostly carried by the acoustic 
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Figure 4: (Color online) Temperature profiles with energy- and momentum-conserving noisy dy- 
namics for a harmonic chain with even number of particles. Other parameters were taken to be 
ma = 0.5, rub = 1.5, Jl = 1R = 1-0 and Tl = 2.0, Tr = 1.0. 

modes. 

IV. SIMULATION RESULTS ON THE EFFECT OF NOISE IN THE DYNAMICS 

As discussed in the introduction, temperature oscillations have been observed in anhar- 
monic chains, where however the oscillations decay with system size. This is expected since 
anharmoncity leads to interactions between phonons which helps to establish local thermal 
equilibrium. A simple model which incorporates phonon-phonon interactions was intro- 
duced in |T3| HI] where the determinsitic dynamics of the Harmonic chain is stochastically 
perturbed. Here we have carried out simulations with this noisy dynamics and looked at 
it's effect on the temperature profiles of the alternating mass chain. There are two cases to 
consider: 

(a) Momentum conserving noise: Here, in addition to the Hamiltonian dynamics without 
pinning, one introduces random exchange of momentum between nearest neighbor particles, 
which occurs with a rate A. This conserves both momentum and energy. In Figs, 
show the effect of momentum conserving noise on the temperature profiles for chains of even 
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we 




Figure 5: Temperature profiles with energy- and momentum-conserving noisy dynamics for a 
liarmonic cliain witli an odd number of particles. Other parameters were taken to be rua = 0.5, 
rub = 1.5, JL = 1R = 1.0 and Tl = 2.0, Tr = 1.0. 




Figure 6: (Color online) Temperature profiles with only energy-conserving noisy dynamics for 
a harmonic chain with even number of particles. Other parameters were taken to be nia = 0.5, 
rub = 1.5, JL = 1R = 1.0 and Tl = 2.0, Tr = 1.0. 
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Figure 7: (Color online) Temperature profiles with only energy-conserving noisy dynamics for a 
harmonic chain with an odd number of particles. Other parameters were taken to be rua = 0.5, 
rub = 1.5, JL = 1R = 1.0 and Tl = 2.0, Tr = 1.0. 

and odd number of particles. In the even case we see that, on introducing noise, the size 
of the oscillations has decreased and the phase of the oscillation on the left half has changed 
sign. For the odd case, the choice of parameters (7/, = 7/? = 1) corresponds to a case with 
no oscillations when A = 0. On introducing noise, A > 0, one gets oscillations very similar to 
the even case. We also see that the oscillation amplitude becomes smaller on increasing 
system size, for both even and odd cases. 

Thus we see that the temperature profile in this system with energy-momentum conserv- 
ing noisy dynamics shows the following qualitative features : (i) The oscillations decay as 
we go into the bulk, (ii) There is a phase shift in the sign of the oscillation amplitude as 
one crosses the center of the chain. The lighter particles at the hot end are always hotter 
than the heavier particles. At the cold end, the heavier particles are hotter. Thus this is 
qualitatively different from the harmonic case, (iii) For large A^, the temperature profile is 
not sensitive to whether A^ is even or odd. These same features have also been observed 
earlier for the alternate mass FPU chain [10] which had a quartic interparticle interaction 
potential (in addition to the harmonic one). The FPU system has momentum conservation 
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and does not satisfy Fourier's Law as is also the case for the system with noisy dynamics. 

(b) Momentum non-conserving case: When the noise only conserves energy but not 
momentum, as can be obtained by randomly reversing the velocity of the i^^ particle at 
rate A, then, as is seen in [15j, the NESS for N oo corresponds to a local equilibrium 
state. This ensures that the Tj in the bulk is the same for i odd or even independently of 
whether is even or odd. In Figs, 
dynamics on the temperature profile for harmonic chains with odd and even number of 
particles. We observe that for the even case, the oscillations in the temperature decreases 
considerably on introducing the noise, and this reduction is greater when is larger. For 
the odd case however, for small systems, introduction of noise produces small oscillations 
in the temperature profile, but these oscillations eventually decrease as the system-size is 
increased. For both even and odd total number of particles, the decay of the oscillation 
amplitude with system size is faster than for the momentum-conserving case and we quickly 
get a linear temperature profile in the bulk of the system. 

V. DISCUSSION 

In this paper we obtained exact expressions for the temperature-profile and the heat 
current in the alternate mass chain connected to heat baths at different temperatures in 
the limit of infinite system size. This proves rigorously that the temperature oscillations of 
successive particles in the bulk persist even in the thermodynamic limit. 

We provided an understanding of these oscillations by noting that in any given normal 
mode, the mean kinetic energy of a particle depends on its mass. In an acoustic mode, the 
heavier particles have higher mean kinetic energy than the lighter ones, while in an optical 
mode, the lighter particles have higher kinetic energy. On connecting the chain to heat 
reservoirs each of the modes are excited to different degrees, depending on the parameters. 
The kinetic energy of a particle gets contributions from all the modes, both acoustic and 
optical and the net result depends on the distribution of energy in the different modes. If 
both the baths have the same temperature, we have an equilibrium steady state in which 
each mode has the same average energy (equipartition). In this case the temperature at all 
sites are equal. The same is true locally when the system is in local equilibrium. 

The situation is different in the non-equilibrium case where we do not have local equilib- 



6ff), we show the effect of addition of velocity flipping 
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Figure 8: (Color online) Simulation results for temperature profile for a two-dimensional N x W 
strip of harmonically coupled particles with a periodic arrangement of masses. The sites on the 
strip are labeled (i, j) with i = 1, . . . ,N and j = 1, . . . , W. Particles at sites with even i + j have 
mass nia and others have masses nn,. Heat baths are attached to all sites on the layers i = 1 
(temperature T^) and i = N (temperature Tr). Periodic boundary conditions are imposed in the 
transverse (j) direction. Upper plot shows the average temperature on succesive layers for chains 
of lengths = 32 and = 64. There are oscillations in the transverse direction also and this is 
shown in the lower plot which shows the temperatures Tij on all sites of a section of the = 64 
chain. Note that from symmetry we have Tij = Tij^2, and this can be observed here. The width 
of the strips were taken to be = 4. The other parameters were taken to be = 0.6, rrib = 1.4, 
7L = 7R = 1.0 and Tl = 2.0, Tr = 1.0. 

rium and there is no equipartition of energy between the different modes. We then expect 
generically that the mean kinetic energy (temperatures) obtained by adding the contribu- 
tions of all modes will depend on the mass of the particle. It is therefore not so surprising 
that we get different kinetic energies for the different masses. From the above explanation we 
expect that temperature oscillations should also occur in higher dimensional periodic har- 
monic systems. Simulation results for two-dimensional strips (see Fig. ^ suggest that this 
is the case, but more extensive studies are necessary to establish the role of dimensionality. 

As already noted there will be no oscillations in the bulk if the NESS is in local thermal 
equilibrium. To achieve this one introduces interactions between the phonons. Interactions 
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between phonons can be introduced for example by adding stochasticity in the dynamics and 
we studied this case numerically. We find that in this case the oscillations are quahtatively 
different from the purely harmonic case and do not survive in the limit of large system size. 
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Appendix A: Details of calculation 



Here we give more details of the derivation for the temperature profile and the current. 
We basically need to evaluate the integral 

^i^L I J, . , 2 ir^ + r. .m2 
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dujoo' IGjHr , (Al) 

where = [— Mw^ + $ — tuF]^^. We consider the case with k = 1. Let us define Ai^m as 
the determinant of the sub-matrix of [— Mo;^ + $ — iuF] that starts from the l^^ row and 
column and ends in the m}^ row and column. We also define Di ^^ as the determinant of the 
sub-matrix of [— MQ^ + $] starting from the l^^ row and column and ending in the m}^ row 
and column. In terms of these one has: 



G+(f]) 



CO 



with Ai_;_i = -Di,i-i - iu;-fLD2,i-i 



(A2) 



Al+l,N — Di+i^N — i^lRDl-^i^N-l 

Ai,Ar = -Di,jv - iw{'^rDi^n-i + 1lD2,n) - uj'^JlJrD2,n-i ■ (A3) 
Let us now define f{l) — Di^2i and g{l) — Di^2i-i- These satisfy the recursion relation. 



"/(/)" 


= B 


"/(/-I)" 









where B 



(2 - m„a;2)(2 - rribOj'^) - 1 -(2 - nibOJ^) 
(2 - m„a;2) -1 

with the initial condition /(O) = 1 and ^'(0) = 0. Hence we get 



"/(/)" 


= B' 


1 


_9{l)_ 








(A4) 
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The matrix B has unit determinant and can be expressed in terms of the Pauh spin matrices 
a as follows: 



B = cosg 1 + la.nsmq = e**^-"-^ 

, ^B {2 - mauj^){2 - TTihUJ^) - 2 

where coso = ir — = 

^2 2 

and r? is a three dimensional unit vector. Hence we get: 



(A5) 



B' 



^^a.n Iq ^ ^^^^^^^ ^ ^ ^-^^^^^ 



B — cos q 1 
sing 



(A6) 



Combining Eq. (A6) and Eq. (A4) we have : 



sin(/ + l/2)q 
sin(g/2) 



sin(/g) 2n 
g{l) = —^{2-maUj') 
smq 



(A7) 
(A8) 



Note that for odd-dimensional matrices with the first mass equal to m^,, the determinant 
would be given by Eq. (A8 ) with replaced by m^. Using these expressions in Eqs. ( A2|A3 ), 
we then get the following forms for the integrals Jj, depending on whether N is even or odd. 
Case(l) - even N: 



'■odd i 



In 
where 



TT 



2 {N-i + l)q 



i2-m,coy + jlu'' 



^2 (N-i)q 

sm2(g/2) 



duuj 



■ 2 (N-i + l)q 



Al,7v| 

:„2 (.N-i)q 



2 sin^(g/2) 



^{2-maOO^) 



2\2 



sin q 



|Ai 



N\ 



Ai, 



N\ 



sm 



(A^+i)g 



sin(g/2) 



sin(g/2) 



+ too 



(A9) 
sin(A^g/2) 



smg 
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Case(2) - odd A^: 



In^d i = / duu - 



TC 



2 (iV~^ + l)g 



2 (iV-Qg 



lAi -12 



TT 



2 (iV-i + l)g . 2 (iV-i)? 

Sm ^ — ^ , 2 2 ^i'^ ^ — 2^ 



/r> 2N2»1"^^ T ^ I 2 



sin^(g/2) 



where 



TT 



^^lA 12 

l^l,Af| 



• (N+l)q 

(2 - m^w^) ^ 7L7ijw^(2 - mbU^) 

smg 



. (7V-l)g 



sm q 



+ 2W ^7l + 7r j 



(AlO) 
sin(iVg/2) 



smg 



We now consider points in the bulk such that x = i/N and (A^ — i)/N remain finite in 



the A^ — > 00 hmit. We now note that, for real values of < g < vr, Eq. (A5) has two allowed 
solutions for u, namely: 

1 



rriamb 
rriamb 



-[l-0(g)] 

1 + <P{q)] 



where 0(g) = [1 — 2mamf,(l — cos g)]^^^ , , < g < vr . 

These correspond to the frequencies in the acoustic and optical branches of the lattice with 
the frequency ranges < w_ < a/2/M and < uj+ < ^/2/{mM), where m (M) 



is the smaller (larger) of the two masses. For frequencies outside these ranges, Eq. (A5) 
gives imaginary values of g. This means that, for these frequencies, terms such as sinA^xg 
grow exponentially with A^. Hence it is clear that, in the limit A^ — i- 00, the integrals in 
Eqs. ( A9|A10 ) only get contributionsfrom frequencies in the acoustic and optical bands. 
Thus for each of the integrals above, we get: 

/2/Af f-^flfirnM) 



dujF{uj) 







dq\ 
dq 



dco- F{uj 
du 



'2/m 



- Fiu^q)) 



duj+ F{uj- 
doj^ 



dq 

dq 



Fiu4q)) 



We now note from Eqs. ( A9|A10 ) that the required integrands F(uj) have factors of the 
form sin^(A^a:;g) in the numerators and Ai jv in the denominators. In the limit N ^ 00 
the factors sin^(A^a;g) in the numerators can be replaced by 1/2. Next we note that the 
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determinant Ai at always has the following form: 



Ai,Ar = A{q) sin(A^g) + B{q) cos{Nq) 



(All) 



where A and B are smooth complex-valued functions. We now obtain the following result 
for any function g{6, (p) which is periodic in both variables: 

r.27r(7V/2) , 



N^oo 



lim / de g{e, NO) = lim 



d(j)g{ 



lim 



i=(N/2) ^2ni 



1 



i={N/2) 



i=l J2iT{i~l) 



dcl>g{- ct>) = - J2 



i=l 



N 



i={N/2) ^ 2n 

lim > — / dipg( 



.27r(i - 1) 



N 



2tt 



2tt 



de / dijg{9,ip) . 



Using this we obtain: 



Ciq) 



\A{q) sm{Nq) + B{q) cos{Nq)\^ 



C{q)dq 



dq- 



27r 



2tt 



dip 
C{q) 



1 



\A{q) sin-?/^ + B{q) cos^^P 



\A{q)B*{q) - A*{q)B{q)\ ^^^'^^ 

Using this we get the asymptotic forms of the various integrals in Eqs. ( A9[A10 ), and these 
lead to the results given in Eqs. ( TTp2 ). When 7l = 7_r = 7, we can explicitly carry out the 
integrals appearing in these expressions and we get the following results. 
Case (1) - even N: 



2(1 + /3) + + 25)2/3 + 2^5^(^ + 1)2 ^ 

J-odd — — ^=^= — . 4 



+ 



2 

rrib 



VWTl^/l + 2(3 + 2f3^fi 



1 + 2/3/i 



1 + 2/3 + 2/32/i 
1 + 2/3 



- 1 



|5|(5 + 1)^ 



+ m„ 



2/3/i 



2 /i(l + 2/3/i)yrT52 ''1 + 2/3/iV VTTP 
-(l + /3)+/3(25-52)- 5^(1 -5)2^3^ 



(A13) 



rual 1 + 



v/rT2^Vl + 2/3 + 2/32;x 



+ /3 



\s\{i-sy 



v/rT52(l + 2/3/i) 



+ 



\6\ _^ yi + 2^ + 2^2^ 



(A14) 



l + 2^fiy VT+^ ^fY+2^ 
= AT^^-^^^^[2/3 + 1 + 2/32;x(l + 52 - |5|v^rT5^) - ^(2/3 + 1) (2/3 + 1 + 2/32;,)] 

(A15) 



where fi = 2mamb , S = rria — mi, , (3 = 72/ {rriamh) . 
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Case (2) - odd N: 



where 



^odd — leven — i (Al6) 



F = ^(B' - 4AC)"\ G = Ct,, H=^-C{l-t,)-A 

rrib rria nib "^a "^a 
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